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Abstract 

We begin an exploration of the physics associated with the general CP-conserving MSSM 
with Minimal Flavor Violation, the pMSSM. The 19 soft SUSY breaking parameters in this 
scenario are chosen so as to satisfy all existing experimental and theoretical constraints assuming 
that the WIMP is a conventional thermal relic, i.e., the lightest neutralino. We scan this 
parameter space twice using both flat and log priors for the soft SUSY breaking mass parameters 
and compare the results which yield similar conclusions. Detailed constraints from both LEP and 
the Tevatron searches play a particularly important role in obtaining our final model samples. 
We find that the pMSSM leads to a much broader set of predictions for the properties of the 
SUSY partners as well as for a number of experimental observables than those found in any of 
the conventional SUSY breaking scenarios such as mSUGRA. This set of models can easily lead 
to atypical expectations for SUSY signals at the LHC. 
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1 Introduction 



The LHC will soon begin operations in earnest and will thereafter begin a detailed study of the 
Terascale. Amongst the many prospective new physics signals that may be realized at the LHC, the 
possible observation of an excess of events with missing energy is awaited with great anticipation. 
Given the strong collection of evidence for the existence of dark matter [IJ and the naturalness of the 
WIMP paradigm [2j, it is no surprise that many beyond the Standard Model (SM) scenarios predict 
the appearance of missing energy signatures at the LHC [3j. Perhaps the most well explored of 
these scenarios is R-Parity conserving, N=l Supersymmetry (SUSY) where the WIMP is identified 
as the lightest super symmetric particle (LSP), which most commonly corresponds to the lightest 
neutralino or gravitino. The simplest manifestation of this picture is the Minimal Supersymmetric 
Standard Model (MSSM), where the gauge group of the SM is maintained and the matter content 
of the SM is augmented only by an additional Higgs doublet. 

While the MSSM has many nice features, perhaps its least understood aspect is the nature of 
the mechanism which breaks the Supersymmetry. There is an ever growing list of possible candidate 
scenarios, including mSUGRA [3j, GMSB [5j, AMSB [6], and gaugino mediated supersymmetry 
breaking [7]. Each of these mechanisms strongly influences both the mass spectra and decay 
patterns of the various sparticles and hence it is highly non-trivial to make model-independent 
statements about the properties of the LSP or about SUSY signatures at colliders such as the 
LHC. This puts forth the question of whether there is some way to study the MSSM in a broad 
fashion while simultaneously adopting as few simplifying assumptions as possible. The purpose of 
this paper is to address this issue. In trying to answer this question, one is immediately faced with 
the daunting task of dealing with the over 120 parameters, arising mainly in the soft-breaking sector, 
which describe SUSY breaking in the MSSM [8j. Clearly the phenomenological study of any model 
with such a large number of parameters is totally hopeless. Obviously some assumptions need to 
be made in order to reduce this large number of independent parameters to a viable subset and 
there are many possible avenues one could follow. Most studies assume a specific SUSY breaking 
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mechanism, employing theoretical assumptions at the GUT scale, and thus reduce the number of 
independent model parameters to 4 or 5. In this investigation, we will not make any assumptions 
at the high scale. We will, however, restrict ourselves to the CP-conserving MSSM (i.e., no new 
phases) with minimal flavor violation (MFV)[9j. Moreover, to simplify matters further and help 
soften the impact of the experimental constraints arising from the flavor sector, we will require 
that the flrst two generations of sfermions be degenerate. This leaves us with 19 independent, real, 
weak-scale, SUSY Lagrangian parameters to consider. These include the gaugino masses Mi^2,3, 
the Higgsino mixing parameter /i, the ratio of the Higgs vevs tan /?, the mass of the pseudoscalar 
Higgs boson m^, and the 10 squared masses of the sfermions (5 for the assumed degenerate first 
two generations and a separate 5 for the third generation). Finally, due to the small Yukawa 
couplings for the first two generations, independent A-terms are only phenomenologically relevant 
for the third generation and will be included in our study for the 6, t and r sectors. This set of 
19 parameters leaves us with what has been called the phenomenological MSSM (pMSSM) in the 
literature [To] and is the scenario that we will study in some detail below. We note that for this 
study we will make the further assumption that the LSP is the neutralino and is a conventional 
thermal relic within the standard cosmology with a radiation dominated era. 

The primary goal of the first phase of this study is to obtain a large set of MSSM parameter 
points that can be used for future studies of SUSY signatures, e.g., signals at the LHC or in 
cosmic rays. To accomplish this, we will perform a scan over this 19-dimensional parameter space 
and subject the resulting specific sets of parameter choices (which we denote as models) to a 
large and well-known (and some perhaps, not so well-known) set of experimental and theoretical 
constraints [11 j. Furthermore, we will perform two independent scans which will employ different 
priors for choosing the above mass parameters; this is important since some bias may be introduced 
into the analysis by how one choose parameter space points. By choosing more than one prior and 
comparing results this sensitivity can be estimated. We will show that while both scans yield 
qualitatively similar results, quantitative differences will be observed in some cases. 
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We will perform our comparison to the experimental data with some care to identify sets 
of parameter space points which are at present phenomenologically viable. As we will see below, 
applying these constraints can be complicated and is not always as straightforward as employed in 
some other analyses. This is particularly true in the case where experimental bounds are predicated 
on specific SUSY breaking mechanisms. As we are not a priori interested in any particular region 
of parameter space {e.g., we are enforcing the WMAP dark matter density[12j observations only as 
an upper bound) or in obtaining a set of points that provide a 'best fit' to the data, we will not need 

to make use of Markov-chain Monte Carlo techniques [HI [Hjliil We will then perform a detailed 
study of the properties of these surviving models and, in particular, highlight those aspects which 
are not generally encountered when only a given set of SUSY breaking scenarios is considered. This 
allows us to get a perspective on the wide range of predictions for various observables that can arise 
in the more general pMSSM framework in comparison to those found in specific SUSY breaking 
scenarios. 

In the next Section, we will discuss our scanning techniques, the ranges of the soft breaking 
parameters that we employ, as well as the values of the SM input parameters that we use. In 
Section 3, we will detail the various theoretical and experimental constraints at some length and 
will explain exactly how they are employed in our analysis. Section 4 presents our analysis and a 
survey of the results, while Section 5 contains an overview discussion of the lessons learned from 
this analysis. Our conclusions can be found in Section 6. 

2 Scans, Parameters Values and Ranges 

To perform our multi-dimensional parameter scan, we need to choose suitable ranges for the soft- 
breaking parameters that we employ and determine how the values of these parameters are picked 

"The MCMC approach asks a different question than we do here, namely 'What regions of the parameter space 
provide the best fits to the data?', by using, for example, a chi-squared analysis. This requires one to assume 
precise knowledge of the probability distributions for both experimental results as well as theoretical predictions. For 
example, given the uncertainties associated with the g-2 of the muon, is it justified to treat the error associated with 
the apparent difference with the SM prediction as Gaussian? These types of questions introduce prejudice into the 
analysis, which we try to avoid, and so we have instead adopted the scanning approach. 
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within these assumed ranges. In the analysis presented here, we will make two independent pa- 
rameter scans. In the first case, we randomly generate 10^ sets of parameters, assmning flat priors, 
i.e., we assume that the parameter values are chosen uniformly throughout their allowed ranges. 
With an eye toward the experimental constraints to be discussed in the next Section, for the flat 
prior scan we employ the following ranges: 

100 GeV <mj<l TeV , 

50GeV < |Mi,2,/u| < ITeV, 
100 GeV < Ms < 1 TeV , 

\Ab,t,r\ < ITeY , (1) 
1 < tan /? < 50 , 
43.5 GeV < < 1 TeV . 

Note that absolute value signs are present in these quoted ranges as we will allow the soft-breaking 
parameters to have arbitrary sign. Given that the soft parameters Mi,Ai and /i can all appear 
with arbitrary signs (i.e., phases) but that only six relative signs are physical[l5] we choose M3 to 
be positive. This scan will generate SUSY spartners with light to moderate masses. 

For our second scan, we randomly generate 2 x 10^ sets of parameters assuming log priors 
for the mass parameters. We then modify the allowed soft parameter ranges as follows: 

100 GeV < mj < 3 TeV, 

10 GeV < |Mi,2,At| < 3 TeV, 
100 GeV < Ms < 3 TeV, 
lOGeV < |A,t,r| < 3TeV, (2) 
1 < tan /? < 60 , 
43.5 GeV <mA< 3 TeV. 

Here, / = Ql, Qs, Li, L3, ui, di, tts, ds, ei, and 63. This expanded range of parameters will allow us 
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access to both very light as well as some heavy sparticle states that may only be produced at the 
SLHC. The goal of this second scan is to make contrasts and comparisons to the flat prior study in 
order to determine the dependence of the model properties on the scan assumptions. The possible 
bias introduced by choosing a particular prior in analyses such as ours is an important problem 
that has been often discussed in the literature [22j. The only way it can be properly addressed is 
to take more than one assumed prior, compare the results and look for differences. For the case of 
our two priors we will see that they yield qualitatively similar results, but the detailed predictions 
in the two cases will be found to be quantitatively different in several aspects. 

Note that the above mass parameters are assumed to be evaluated at the SUSY scale, 
i.e., by convention the geometric mean of the two stop masses ^ 1 TeV. The physical spectra for 
the sparticles themselves are generated in all cases using the code SuSpect2.34|10j. Once a set of 
parameters is chosen, we subject the resulting MSSM model to the entire set of constraints discussed 
in the next Section; those models which survive will be analyzed in further detail below. Of course, 
as models are passed through each successive experimental 'filter,' the number of surviving models 
will be further and further reduced. 

Another important aspect of such a broad study is the set of input values for the SM 
parameters, e.g., as{Mz), mt, M\y, etc., that we employ. As is well-known, some of the predictions 
for the various observables we analyze below can be sensitive to these particular numerical choices. 
While we have chosen these numerical values with some care, we are not advocating any one 
particular set of values for these parameters. We provide their values here only for completeness 
since it is important for the interested reader to know exactly how our analysis was performed. 
Table [1] lists the values of these SM parameters used in our analysis. While we only quote central 
values here, many of the relevant errors are included when evaluating the set of observables we 
discuss in the next Section. We follow this approach, rather than including the SM input values in 
our scan, since they are essentially known quantities unlike the pMSSM parameters. 
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Parameter 


Value 


aiMz) 


127.918 [Ref.[16j] 


as{Mz) 


0.1198 [Ref.tlTj] 


Mz 


91.1875 GeV [Ref.[IS]] 


Tz 


2.4952 GeV [Ref.[l6l US]] 


sin ^-ii; Ion— shell 


0.22264 [Ref.p]] 


Mw 


80.398 GeV [Ref.[l8]] 


Tw 


2.140 GeV [Ref.[l8]] 


GeV) 


128 MeV [Ref.[l6]] 


pole 


1666 MeV [Ref.ilTj] 


mb 


4.164 GeV [Ref.flT]] 


pole 


4.80 GeV [Ref.dZ]] 


pole 


172.6 GeV [Ref.[l9]] 




0.2255 [Ref.[l6]] 


Vch 


41.6 X 10-3[Ref.pJ] 


Vub 


4.31 X 10-3[Ref.[16j] 




0.104 [Ref.[l6]] 


ruBa 


5.279 GeV [Ref.flB]] 




216 MeV [Ref.[l6]] 




1.643 ps [Ref.[l6]] 




230 MeV [Ref.[20]] 


TBs 


1.47 ps [Ref.[2l]] 



Table 1: Central values of the SM input parameters used in our analysis. 
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3 Theoretical and Experimental Constraints 



Apart from the restrictions imposed by the ranges chosen above for the 19 pMSSM soft breaking 
parameters, there are a large number of both theoretical and experimental constraints that need 
to be considered in order to obtain phenomenologically viable models. Many of these are familiar 
while others will be somewhat less so to the arbitrary reader. The direct application of subsets of 
these numerous constraints and the specific ranges which are considered viable varies widely in the 
literature from one analysis to another [22]; hence, we will discuss them here in some in detail. 

The first set of constraints we impose are theoretical and rather standard and are applied 
while generating the sparticle spectrum with the SuSpect code: the sparticle spectrum must be 
tachyon free and cannot lead to color or charge breaking minima in the scalar potential. We also 
require that electroweak symmetry breaking be consistent and that the Higgs potential be bounded 
from below. Note, however, since we do not impose grand unification (as we are ignorant of the 
ultraviolet completion of the MSSM if we do not implement a specific SUSY breaking scenario), 
we do not require that the RGE evolution of the couplings be well-behaved up to the GUT scale. 
Because we are only concerned with the physics between the weak and TeV scales and the overall 
effects of RGE evolution in this range will be rather modest, it is unlikely that a Landau pole or 
other irregularities can be generated in this narrow mass window. 

At this point in our analysis, we will make the further assumption that the WIMP be a 
conventional thermal relic so that the LSP can be identified as the lightest neutralino; the possibility 
of the LSP being the gravitino or the axino will be ignored in the rest of the present analysis [23 j. 
The scenario where the gravitino is the LSP [2^ will be discussed in a future analysis within a 
different astrophysical framework. We note that the case where the sneutrino is the LSP can be 
easily eliminated in the pMSSM by combining several of the constraints discussed below. The 
possibility that the LSP is a non-thermal relic can also be entertained and may lead to somewhat 
different results. 
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The code micrOMEGAs2.21 |25j takes the MSSM spectrum output from SuSpect and ahows 
us to implement the constraints arising from a number of precision and flavor measurements: 
precision electroweak constraints via Ap, the rare decays b ^ and Bs — > as well as the 

g — 2 oi the muon. We take the 95% CL allowed experimental range for Ap directly from the 
analysis presented in the 2008 PDG[16]: —0.0007 < Ap < 0.0026. For 6 — > 57 we use the combined 
experimental result from HFAG[21J, Bb^s'y = (3.52±0.25)-10^^, as well as the recent SM theoretical 
predictions as given by both Misiak et al.[26\ and Becher and Neubert[27j. Combining both the 
experimental and theoretical errors, we require, rather conservatively, that the predicted branching 
fraction lie in the range Bb^s'y = (2.5 - 4.1) • 10"^. For the decay B^ p^p , we employ 
the recently reported combined limit on the branching fraction obtained by CDF and D0|28j. 
BBs^fMf, < 4.5 • 10^8 at 95% CL. Note that our philosophy will be to apply some experimental 
constraints rather loosely, as they can always be tightened in obtaining a final sub-sample. This 
will also allow us to easily study the sensitivity of the various observables to the particular values 
of the soft breaking parameters as well as determine the favored ranges of the various observables 
themselves. 

As is well-known, measurements of the anomalous magnetic moment of the muon, {g — 2)^, 
differ p9] from the recently updated prediction of the SM by more than ~ 3a: this difference is 
found to be A{g - 2)^ = (30.2 ± 8.8) • 10"^° in one recent analysis[30] and (29.5 ± 7.9) • 10-^° in 
another [3T]. with both analyses using low energy e^e~ data as input. It has been argued that this 
deviation might be a signal for new physics beyond the SM, e.g., light Super symmetry. In order 
to explore this possibility more fully within the general pMSSM context and also to allow for the 
scenario that the SM value may essentially end up being correct once the theoretical, experimental 
and input errors are more fully understood and reduced, we will implement the loose constraint 
that (—10 < A{g — 2)^ < 40) • 10^^*^ in our analysis. This will also allow us to observe how the 
MSSM predictions are distributed within this range for the set of surviving models. 
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In addition to these constraints which are essentially built into the micrOMEGAs2.21 code, 
we employ several others that arise from heavy flavor measurements: (i) the branching fraction 
for B — > TV] combining the results from HFAG[2T] and those presented at ICHEP08[32] we will 
require this branching fraction to lie in the range Bb^tu = (55 - 227) • 10"^ Here we make use 
of the theoretical analysis by Isidori and Paradisij33j as well as Erikson, Mamoudi and Stal|34j 
which include the SUSY loop contributions to the bottom quark Yukawa coupling, (ii) With 
the assumption of MFV[9] and degenerate first and second generation sfermions, the resulting 
constraints from meson-antimeson mixing on the squark mass spectra are relatively modest [35], We 
performed a detailed numerical study which demonstrated that these conditions can be generally 
satisfied if the ratio of first/second and third generation squark soft breaking masses (of a given 
charge and helicity) differ from unity by no more than a factor of ^ 5. We will impose this 
restriction on the mass spectra of the models we consider. We also apply similar restrictions in the 
slepton sector. 

It is well-known that restrictions on the MSSM parameter space arise from employing the 
LSP as a form of dark matter. Two major constraints are found in this case: the first arises from 
the 5 year WMAP measurement [12] of the relic density, which we employ only as an upper limit 
on the LSP contribution, i.e., Qh^\LSP ^ 0.1210. Here we simply acknowledge the possibility 
that even within the MSSM and the thermal relic framework, dark matter may have multiple 
components with the LSP being just one possible contributor0 The second constraint arises from 
direct detection searches [36] for dark matter via both spin-independent and spin-dependent LSP 
elastic scattering off nuclei in the laboratory. Over essentially all of the parameter space, the 
spin- independent constraints are presently found to be dominant. We include the cross section 
bounds arising from XENON10[37], CDMS[38], CRESST I[39] and DAMA[l0] data. Given the 
uncertainties from low energy physics {e.g., nuclear form factors) in the determination of these 

cross sections, we allow for a factor of 4 uncertainty in these calculations. This numerical factor 

**We note that within our analysis framework, given the order in which our constraints are applied, the use of the 
WMAP constraint removes ~ 50% of the models which have survived up to that point in the case of flat priors. 
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was obtained by varying the low energy input parameters independently over their allowed ranges 
for several benchmark points. Thus we will allow models to predict cross sections as much as 4 
times larger than the usually quoted experimental bounds. We employ micrOMEGAs2.21 |25j in 
evaluating both Qh'^ as well as the dark matter scattering cross sections. We note that a 'low' 
prediction for the thermal relic density in some of our models will lead to an effective direct search 
cross section which is smaller than what one would obtain if the WMAP bound were saturated. 

We now turn to the constraints imposed by LEP data. First, we require that there be no 
new charged sparticles or charged Higgs bosons with masses below Mz/2 due to the lack of evidence 
for such states in both direct or indirect measurements at the Z pole; we also impose this same 
constraint on the lightest neutral Higgs boson. Based on data from LEPIIjJl], we further require 
that there be no new stable charged particles of any kind with masses below 100 GeV. We note 
that such a situation may occur if, for example, the nLSP is charged and highly degenerate with 
the LSP, which is a common occurrence in the AMSB scenario. Additionally, we constrain any new 
contributions to the invisible width of the Z boson to be < 2 MeV|18]. This will restrict the masses 
and couplings of light neutralinos that have significant Higgsino content and thus may appear in 
the decay Z XiXi- While it is unlikely that decays such as Z — > X2X1 will be kinematically 
accessible, in such cases the X2 '^i^^ more than likely decay into visible final states. We isolate 
such parameter space model points in the analysis below and scrutinize them furtheiP^ (using the 
SUSY-HIT1.1[42J to analyze the X2 decays). Furthermore, we note that for the range of sfermion 
mass soft breaking parameters we consider, Z decay to pairs of sneutrinos is not kinematically 
allowed so that this final state cannot contribute in any way to the purely invisible width. 

Data from LEPII provides further direct search constraints on the sparticles, although 
they are oftentimes sensitive to detailed assumptions about the model; we provide only an outline 
of these many bounds here. ALEPH[43j has placed a lower limit on the light squark masses, 

assuming that the gluino is more massive than the squarks, via their decay to a jet+LSP, i.e., 

^^In the one model that falls into this category we find that X2 essentially stable but its contribution to the 
invisible width of the Z is tiny. 
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jet + missing energy. Provided that the mass difference between the squark and the LSP, Am, 
is > 10 GeV to avoid very soft jets, they obtain a lower hmit of 92 GeV on the squark masses. 
Given our soft breaking parameter spectrum ranges above, we employ this constraint directly, 
including the Am cut. Models in our scan with Am < 10 GeV are hence not affected by this search 
limit. A comparable search [31] for light sbottoms using the same sort of decay pattern results in 
a lower bound of 95 GeV on the mass of this sparticle. Corresponding lower bounds have been 
placed [35] on the masses of right-handed sleptons undergoing decays to leptons plus missing energy 
of m ^ 100(95,90) GeV for the selectron(smuon,stau). This limit is only applicable if the slepton 
masses are at least a few percent larger than that of the LSP, otherwise the final state leptons will 
be too soft. As in the case of squark production, this mass splitting requirement is critical and is 
included in our analysis. Strictly speaking, these bounds are only applicable to the right-handed 
sleptons, however, they may also be applied to left-handed sleptons provided the corresponding 
sneutrino t— channel exchange contribution is not very important. We make this assumption in 
our analysis. An analogous situation applies to chargino production. If the LSP-chargino mass 
splitting is Am > 2 GeV, a direct lower limit of 103 GeV on the chargino mass is obtained from 
LEPII data. However, if this splitting is Am < 2 GeV, the bound degrades to 95 GeV, provided 
that also Am > 50 MeV, otherwise the chargino would appear as a stable particle in the detector 
and would then be excluded by the stable particle searches discussed above. In the case where the 
lightest chargino is dominantly Wino, this limit is found to be applicable only when the electron 
sneutrino is more massive than 160 GeV. There are also two searches for light stops to consider: if 
ii is too light to decay into Wbxi, then the search for the decay ii — > cxi provides a lower bound 
on the stop mass of 97 GeV. If ti is heavy enough to decay into £bi), then the corresponding lower 
bound is 95 GeV. 

As provided in detail by the LEP Higgs Working Group [37], there are five sets of constraints 
on the MSSM Higgs sector imposed by LEPII data. We employ the SUSY-HIT routine to analyze 
these, recalling that the uncertainty on the calculated mass of the lightest Higgs boson is approx- 
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imately 3 GeV[l6] as determined by SuSpect. The first pair of these constraints apphes to the 
products g\zh ^ -^(^ ~^ bb,T^T^), where gzzh represents the ZZh coupling, relative to their SM 
values as a function of the light Higgs mass. A second pair of constraints applies to the correspond- 
ing couplings QzhA ^ B{hA bbbb, bbr^T^) as a function of the sum of the h and A masses. The 
last of these LEPII Higgs bounds applies to the mass of the charged Higgs boson as a function of its 
branching fraction into the tv and cs final states. We include all of these experimental constraints 
in detail in our results presented below. 

The final set of constraints we incorporate arises from Tevatron data. We first consider the 
restrictions imposed on the squark and gluino sectors arising from the null result of the multijet 
plus missing energy search performed by D0|48j. In our study, we generalize their analysis to render 
it model independent. For each of our models that have survived the set of constraints thus far, we 
compute the NLO SUSY cross sections for squark and gluino production using PROSPINO2.0[l9]. 
Once produced, these sparticles are decayed via SDECAY/HDECAY to obtain the relevant decay 
chains and branching fractions and are then passed to PYTHIA6.4[50j to include hadronization 
and fragmentation. We then use PGS4[5T] to simulate the DO detector and to impose the analysis 
kinematic cuts; we tune PGS4 by reproducing the results and efficiencies for the three benchmark 
mSUGRA points chosen by DO in their published multijet study. Given an integrated luminosity of 
2.1 fb~^, we find that the 95% CL upper limit on the number of signal events from combining all 
of the production channels is 8.34 using the method of Feldman and Cousins[S2]. If our calculation 
for a model leads to an event rate larger than this value, that model is removed from the remaining 
sample. 

Analogously, we employ constraints from the CDF search for trileptons plus missing energy [53] 
which we also generalize to the full pMSSM. Our procedure is essentially the same as in the jets 
plus plus missing energy analysis described above, except that we now employ a CDF tuning for 
PGS4 which we obtain by reproducing the CDF benchmark point results. Here we use the leading 
order cross section together with a universal K-factor of 1.3 to mimic the full NLO cross section. 
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We used PROSPINO to check that this K-factor is a good approximation for the range of gaugino 
masses relevant to the Tevatron. In this case we only make use of the '3 tight lepton' analysis from 
CDF as it is the cleanest and easiest to implement with PGS4. The 95% CL upper bound on the 
possible SUSY signal in this channel is then 4.65 events assuming a luminosity of 2.02 fb^^ as used 
in the CDF analysis. Again, parameter sets leading to larger event rates are dropped from the 
remainder of our analysis. 

Over the narrow mass range 90 < niA < 100 GeV, searches for SUSY Higgs signals by 
both CDF and D0[5l] provide the approximate constraint tan/3 > 1.2mA — 70. We include this 
restriction in our analysis. 

Finally, both CDF^55j and D0^56j have placed lower limits on the masses of heavy stable 
charged particles. In our analysis we employ the stronger DO constraint which we take to be 
of the form > 206|C/i^p + 171|C/ih,p GeV at 95% CL in the case of charginos. Here, the 

matrix U determines the Wino/Higgsino content of the lightest chargino and is used to interpolate 
between the separate Wino and Higgsino results as quoted by DO. As we will see below, this is an 
extremely powerful constraint for a large number of our models since we find that chargino-LSP 
mass degeneracies are common in our model sample, particularly when the LSP is nearly pure 
Wino or Higgsino or a combination of these two cases. 

We note in passing that both CDF and DO also have a set of analyses aimed at searches for 
light stops and bottoms [57j with various assumptions about the SUSY mass spectrum and possible 
decay channels. Most of these searches are only applicable if these sparticles are lighter than the 
top quark. Since these are a difficult set of analyses to generalize to the pMSSM we will constrain 
our model set to the subset of those which contain stops and sbottoms with masses larger than 
TTif. We separately store those models satisfying all of our other constraints above but contain 
lighter stops and sbottoms (with masses above the LEPII limits) for a further later detailed study. 
However, we will find that set of models in this category are found to be extremely rare as will be 
obvious from the stop and sbottom mass distributions presented in the next Section. 
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A final comment on the Tevatron SUSY searches is in order. As can be seen above, there are 
only two generic, flavor- independent, high Et searches for SUSY particles at the Tevatron. Given 
our large parameter space we will see that these searches are relatively easy to evade. Some fraction 
of the remaining parameter space may, however, be closed if further searches of these kind can be 
performed at the Tevatron. An obvious example of this is the usual 'cascade' decay picture which is 
a hallmark of the LHC SUSY searches where, e.g., squarks decay via electroweak gauginos leading 
to a multi-jet, multi-lepton plus missing Et final state. Such a decay signature may certainly occur 
for many of our models and it thus remains only a question of the signal rates vs. backgrounds 
and the development the most efficient cuts to determine how effective such a search might be at 
the Tevatron. We encourage our Tevatron colleagues to consider performing such a search and to 
report their findings in as model independent a manner as possible. 

4 Analysis and Results 

We now describe the results from our analysis for both the cases of flat and logarithmic priors. 
4.1 Flat Priors 

In the case of flat priors, we randomly generated 10'' parameter space points (i.e., models). After 
subjecting these models to all of the constraints detailed above, we found that only ~ 68.5 ■ 10^ 
satisfy all of the restrictions. Here, we discuss the properties of these surviving models in some 
detail. The first characteristic one would like to know is the nature of the sparticle spectra in these 
various models. Since we cannot analyze so many complete spectra here, we will instead examine 
a number of specific properties of the models which then will provide an overall qualitative feel for 
their general behavior. 

Figures [T]- [6] show histograms of the distributions of the various sparticle masses, as well as 
for the charged and neutral Higgs boson masses, in our set of surviving models . In Fig.[Tl we display 
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the distributions of the selectron/smuon, stau and sneutrino masses. We observe several features: 
(i) right-handed selectrons/smuons tend to be Ughter than their corresponding left-handed sleptons, 
{ii) the lighter stau tends to be the lightest charged slepton, i.e., there is an inclination for it to be 
lighter than selectrons or smuons with either helicity label; this is most likely due to stau mixing 
effects, (Hi) tau sneutrinos have a proclivity to be lighter than electron/muon sneutrinos (which is 
not surprising, recalling the connection between the masses of the left-handed charged sleptons and 
their associated sneutrinos), and (iv) a reasonable fraction of these models, ^ 10%, would predict 
that some type of slepton is kinematically accessible at a 500 GeV e~^e~ linear collider. 

We present the squark mass distributions in Fig. [2l where we see that the lightest stops and 
sbottoms are naturally heavier than the top mass (and thus the constraint on stop/sbottom masses 
imposed above is found to have little effect). The upper left panel reveals particularly interesting 
results for the first/second generation (i.e., 'light') squarks. Here we see that there are a significant 
number of models with first/second generation squarks masses below ~ 300 GeV, i.e., well below 
the limits usually quoted from Tevatron squark searches[48]. It is important to remember that all 
the models shown here have been subjected to and have passed the Tevatron search constraints from 
DO. However, the Tevatron searches are targeted at mSUGRA models, which have a very restricted 
type of mass spectra. There are several ways that light squarks may survive the DO analysis and 
not be excluded: (i) if a left-handed squark is light it will tend not to decay immediately to missing 
energy but instead undergo a cascade (assuming the electroweak gauginos are also light). In this 
case the DO search strategy can easily fail as the search criterion are invalid, (ii) If a right-handed 
squark is relatively light, it will more than likely decay into a jet plus the LSP. However if the mass 
difference between the squark and the LSP is relatively small (something that does not normally 
happen in mSUGRA) the resulting events will contain rather soft jets and will not be affected by 
the kinematic cuts imposed by DO. (iii) Depending on the gluino mass, the production rate of 
the SUSY signal can be suppressed so that there are not enough events to reach the statistical 
boundary for an observable signal. This result is important as relatively light squarks could exist 
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and may have been missed by the Tevatron searches. In fact, there appears to be a reasonable 
window for hght right-handed squarks to be kinematicahy accessible at a 500 GeV linear collider. 
Hunters of squarks at the LHC must keep this possibility in mind when constructing detailed search 
strategies. 
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Figure 1: Distribution of slepton masses for the set of flat prior models satisfying all of our con- 
straints: selectrons/smuons(top left), staus(top right) and sneutrinos (bottom). 

Figure [3] displays the mass distributions for both neutralinos and charginos. While the LSP 
tends to be quite light, we see it is also common for both and X2 t° ^^^^ t° quite light with 
a mass almost always below 250 — 300 GeV. As we will see below, the almost pure Higgsino 
LSP scenario is reasonably common in our model set and the rough degeneracy of these three 
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states then follows nearly automatically. The almost pure Wino LSP scenario is also quite common 
amongst our models and it also predicts a near LSP-x^^ degeneracy. Both of these scenarios can 
lead to long-lived charged particle signatures at the LHC. The mass distribution for gluinos, and 
its correlation with the LSP mass, is given in Fig. HI Here, again, we see that many models allow 
for a rather light gluino, something we would not have expected based on mSUGRA due to the 
results of the DO jet plus missing energy search. However, most of these light gluinos escape the 
DO search due to, e.g., the smaller mass splitting with the LSP as discussed above in the case of 
squarks; this leads to softer jets that avoid the DO search requirements [58]. 

In Fig. [5] we see that the A, H and masses tend to track each other even for rather 
low mass values. A few percent of our models have heavy Higgs bosons which are reasonably light 
and are thus in the non-decoupling scenario. We note that there is a small depletion in the heavy 
Higgs counts when Tnjj ~ 2mh. This is due to a small problem with the calculation of the 
H ^ 2h decay width near threshold [62] in the version of SUSY-HIT that we employ. In order to 
avoid this problem, we have removed a small region of the phase space, leading to the depletion 
observed here in the H, A and distributions. Figure [6] shows the distribution for the lightest 
Higgs mass (recall there is a 3 GeV uncertainty in the SuSpect calculation). Here we see that (i) 
it is difficult to generate masses much above 125 GeV (this is most likely related to our chosen 
parameter range that stop masses be below ~ 1 TeV), and (n) a small set of models predict light 
Higgs bosons with masses as low as 85 — 90 GeV. Clearly in such models the Higgs couplings are 
either quite different than in the SM with, e.g., reduced ZZh couplings with a lighter itla (such as 
in the non-decoupling region), and/or the light Higgs decays into LSP pairs with a large branching 
fraction (while simultaneously avoiding the constraint on Z decay into invisible particles). In fact, 
we find that both situations do occur in our set of models. We note in passing that we encounter 
a second SUSH-HIT issue in the case of very light Higgs in the mass range between 80 and 90 
GeV. Here, an anomalously large partial width for the 77 mode is returned by the code yielding 
an invalid result. These few cases appear in both the flat and log prior samples and are dropped 
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from our final model sample. 




4500 




[GeV] 

Figure 2: Distribution of squark masses for the set of flat prior models satisfying all of our con- 
straints: first/second generation squarks are in the upper left panel, stops in the upper right panel 
and sbottoms are in the lower panel. 



The properties of the nLSP are examined in Fig. \7\ which shows the various possible iden- 
tities of the nLSP in our surviving flat prior model set, as well as the mass splitting between the 
LSP and the nLSP as a function of the LSP mass. These two properties essentially determine the 
collider phenomenology of a model. From this figure, we see that the lightest chargino plays the 
role of the nLSP much of the time but many of the other sparticles can also be the nLSP at the 
level of a few percent. The second most common nLSP is found to be X2 occurring in ~ 6 — 7% of 
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Figure 3: Distribution of neutralino and chargino masses for the set of flat prior models satisfying 
all of our constraints. 
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Figure 4: Distribution of gluino masses and a comparison of the gluino and LSP masses for the set 
of flat prior models satisfying all of our constraints. 
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Figure 5: Distribution of heavy and charged Higgs masses for the set of flat prior models satisfying 
all of our constraints. 
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;ure 6: Distribution of the light Higgs mass on a log scale for the set of fiat prior models satisfying 
of our constraints. 
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the cases; these are mainly pure Higgsmo LSP scenarios. Other identities of the nLSP may lead to 
unusual signatures at the LHC. In the lower panel of this figure we see not only the nLSP identity 
but also the nLSP-LSP mass splitting, Am, as a function of the LSP mass itself. Here we see that 
there is a strong concentration of models with small values of Am in our sample; in hindsight this 
is not too surprising since, as we shall see below, the LSP in our model set is commonly either 
pure Wino or Higgsino or a small admixture of these two cases. The bulk of the mostly empty 
square region which appears on the lower left-hand side of the figure has been removed by the 
stable chargino search at the Tevatron; only non-chargino nLSPs remain in this regiorjtt]. Again, 
we emphasize that the LSP in the models with small values of Am tend much of the time to be close 
to a Wino or Higgsino (or combination thereof) electroweak eigenstate as we will see below. We 
note that a stable heavy charged particle search at the LHC should exclude or discover the models 
with heavier Chargino nLSPs and small values of Am. There are many interesting features worth 
noting in the lower panel of this figure. For example, we see that there are a cluster of models with 
light charginos, i.e., masses less than 100 GeV, but with small ( ^ 1 GeV) values of Am. Such 
particles, though sufficiently short-lived, would only leave soft tracks in a LEP detector and would 
likely to have been missed. We also see a set of models with a very small Am where X2 the nLSP 
but with masses such that Z X2X1 is kinematically forbidden. Again, such final states would 
be nearly impossible to observe and would have been missed at LEP. This figure also reveals cases 
with light squark nLSPs and small values of Am. Furthermore, we observe that models in this 
figure with the lightest LSPs and also have larger values of Am ~ 10 — 300 GeV, can have almost 
any SUSY particle as the nLSP. The small gap in the region 1 ^ Am ^ 30 for LSP masses ^ 100 
GeV may be easily explained. One possibility is that in the many models where we encounter an 
almost pure Wino or Higgsino LSP with a small mass difference with a chargino nLSP, the mass 
of such a chargino must be in excess of ^ 100 GeV due to LEP constraints. Bino LSPs, on the 

otherhand, are not subject to such constraints and can be much lighter though are found to be 

"''"'"Note that sleptons with masses in this range have cross sections which are too small to be excluded by the 
Tevatron search. 
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somewhat less common in our sample.. 

Based on these results there are many interesting possible scenarios that can arise at the 
LHC. Here, we give two examples involving small mass splittings between the nLSP and LSP. In the 
first case we consider, we find that there are a class of models with a squark nLSP with very light 
masses ~ 80 — 120 GeV and mass splittings from the LSP of only a few GeV or less. In the low end 
of this squark mass range, such events would have appeared as a two-photon background at LEP, 
whereas for the full mass range the relative softness of the jets in the final state would imply that 
the squarks would have been missed by the DO multijet search. Such sparticles would be notoriously 
difficult to detect at the LHC; it is perhaps possible that using event samples with additional gluon 
ISR would be useful [58] but would require further study. A second interesting scenario is the case 
where X2 nLSP having a small mass splitting with the LSP. The properties of X2 neutralinos 

are usually determined via the cascade decays of heavier colored sparticles produced with large 
cross sections or when they are produced in association with the lightest chargino via the trilepton 
mode discussed above. The most common case where X2 nLSP in our model sample is for 

both it and the LSP to be dominantly Higgsino-like, further implying that the lightest chargino 
would also be nearby in mass. In this situation, due to the Higgsino-like nature, we would not 
expect significant trilepton production and the appearance of X2 cascade decays would be most 
unusual. If X2 is produced in a cascade decay (by, say, stop pair production followed by the decay 
ti X2)> the X2 decay to the LSP would be controlled by the size of the allowed phase space which 
may be well below 1 GeV in some cases leading to a possible secondary vertex. For some mass 
splittings, decays through a virtual Z or h may dominate leading to an additional pair of (b)jets 
or leptons. In other cases, decays into a photon would be dominant leading to a final state with 
additional non-pointing photons. 

Since the nLSP can be almost any SUSY particle, and the corresponding mass splittings 
can be small in all cases, essentially all of the conventional SUSY long-lived particle scenarios can 
be captured in our model set. For example, long-lived stops or staus (as in GMSB[5]), gluinos (as in 
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Split SUSY[59]) as well as charginos (as in AMSB[6]) all occur in our sample as do other possibilities 
including long-lived selectrons, sneutrinos and sbottoms. Of course, long-lived Xi neutralinos can 
occur in GMSB, where they are the nLSP, whereas here long-lived neutralinos can only be the 
X2- As an aside, it is interesting to note that the set of mass spectra obtained above for the 
various sparticles, as well as the large number of models with small nLSP-LSP mass splittings, are 
qualitatively similar to the corresponding ones in a set of models given in Ref. pO] . These models 
were found to be difficult to distinguish at the LHC and were subsequently further analyzed in 
detail in the ILC context [61j. As we will see below, this similarity will remain valid in the case of 
our models generated via log priors as well. This present analysis thus demonstrates that there is 
nothing special (or qualitatively wrong) about that particular [60J model set. 

Figure [8] shows the detailed Wino/Higgsino/Bino content of the LSP in our set of models. 
Here we see that, as advertised, the LSPs in our sample tend to be either a rather pure electroweak 
eigenstate or a mixture of only two of these states. This occurs quite naturally when one of 
the parameters in the neutralino mass matrix takes on a value which is much smaller than the 
others. LSPs with large weak eigenstate admixtures, that would populate the central regions of 
these figures, are seen to be relatively rare in our model sample. In addition, both Wino-like 
and Higgsino-like LSPs are seen to be somewhat (i.e., weakly) dominant over the Bino-like case 
familiar from mSUGRA although all three of these possibilities are very well represented in our 
model sample. 

Complementary information on the Wino/Higgsino/Bino content of the LSP is provided 

in Figure [9l Here we see several things: (i) light LSPs tend to have a large Bino content with a 

possible moderate amount of admixing with the Higgsino, and (ii) for tiny mass splittings with the 

LSP the wino content of the LSP is quite high, while for large splittings the Bino content apparently 

dominates. The Higgsino content is large in the mass splitting region between these two extreme 
cases. 

Let us now turn to the pMSSM postdictions, i.e., the distribution of values from our model 
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Figure 7: (Top) Identity of the nLSP. (Bottom) nLSP-LSP mass splitting as a function of the LSP 
mass, with the identity of the nLSP as labeled. Both assume flat priors. 
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Figure 8: Wino/Higgsino/Bino content of the LSP in the case of flat priors. 
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Figure 9: Fractional Wino/Higgsino/Bino content of the LSP in the case of flat priors as functions 
of the LSP mass(top) and the nLSP-LSP mass splitting(bottom). Note that every model has three 
entries in both of these figures. 
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set for the experimental quantities used as constraints in obtaining that model set. Figure [10] 
displays these postdictions for several of the experimental observables, as well as for the values of 
tan f3, from our surviving model sample. Note that this particular distribution peaks at values near 
tan/3 ~ 12 and favors smaller values of this quantity. As expected, the distribution for 6{g — 2)^ 
is approximately bimodal, depending upon the sign of the Higgsino mixing parameter fi. Here, 
we see that most models do not predict a large enough shift in this observable to fully 'explain' 
the central value of the observed discrepancy with the SM. We also see that while there is a slight 
preference for somewhat larger values of B{b s'j) the predicted values for B[Bs fi^fi^) lie 
significantly below the current Tevatron bound and hover around the SM prediction for this decay. 
Generally, B{B — > tv) is not found to be a great discriminator amongst our model set. Note 
that the distribution for the predicted value of the relic density Qh"^ is peaked at rather small 
values implying that in most of our models the dark matter must be dominantly non-LSP in origin. 
However we also see that there is a long tail leading to a significant number of models with Qh'^ 
values close to the WMAP upper limit. 

We show a number of relevant predictions for dark matter in our flat prior model sample 
Fig. [TTJ In the top panels we see the prediction for il/i^ as a function of both the LSP mass as 
well as the LSP-nLSP mass splitting. Am. Note that the predicted range for Qh"^ is quite large 
but is found to peak near ~ 0.004, which is about 1/30 of the observed value obtained by WMAP 
as discussed above. These models then would require that there be a substantial amount of dark 
matter from other sources such as axions. The range of predictions presented here for the relic 
density is found to be much larger than those obtained by other analyses (.22! that were restricted 
to specific SUSY breaking scenarios. We also show the scaled spin-independent as well as spin- 
dependent LSP-proton scattering cross sections relevant for direct detection experiments. Note 
that again the predicted ranges that we find for these scattering cross sections are far larger than 
those based on mSUGRA, as obtained by, e.g., Ref.[63]. In fact, we find that the range of values 
we obtain from our pMSSM model set essentially covers the entire theoretical anticipated region 
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Figure 10: Distributions of predictions for several observables as well as tan (3 for our model sample 
subject to the constraints discussed in the text. The blue and green dashed lines show the SM 
predictions as well as the current central values obtained by experiment, respectively. 
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found in this reference for all different types of theories beyond the SM. Thus, this observable can 
not be used as a model discriminator, e.5., to distinguish between SUSY versus Universal Extra 
Dimensions or Little Higgs thermal relics. Clearly these results show that the range of expectations 
from the pMSSM is far larger than those predicted within any specific SUSY breaking scenario. 

The amount of fine-tuning in our model sample is described in Fig. [T2l Here we see that 
there is a clear set of pMSSM models with a rather reasonable degree of fine-tuning, unlike the 
current situation within mSUGRA|22j. Following the analysis presented in Ref. [M], we define the 
quantities 

d log m| 



A{0 



(3) 



where ^ are the set of relevant Higgs Lagrangian parameters within the pMSSM. Then, for example 



A{b) 
A{ml) 
A{ml) 



rrir 



l + !^i±^tan^2/3 
mi, 



1 + 



m 



tan^ 2P, 



cos 2(3 + ^ cos^ (3 



rrir 



- cos 2/3 + ^ sin^ (3 
2 



1 



1 + 



COS 2/3 
1 

COS 2/3 



+ ^^i±^tan2 2/? 
mi 



+ 



rriA + mi 



m , 



tan^ 2(3 



(4) 



where it is has been assumed that tan/3 > 1. Here, m^^ are the effective mass parameters for the 

two Higgs doublets generating masses for the u, d-type quarks, respectively, and h is the parameter 
in the Higgs potential that can be rewritten in terms of tan/3 and m^. The overall fine-tuning A 
is defined by adding these four quantities in quadrature. Values of A far above unity then indicate 
large amounts of fine-tuning. In this figure we see that a reasonable set of our surviving models can 
have very low values of finetuning. This figure also shows that the overall amount of fine-tuning 
indeed peaks at rather low values in our set of models, far smaller than that obtained in the case of 
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Figure 11: Distributions of relevant dark matter predictions from our model sample subject to the 
constraints discussed in the text. All cross sections are in pb. 
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mSUGRA. Note that the fine-tuning percentage given in the bottom panel of this figure is simply 
the value of 1/A expressed as a percent. 
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Figure 12: The distribution of fine-tuning (= A) in our set of models and as a function of the LSP 
mass. Also shown is the correlation of the amount of fine-tuning with the predicted dark matter 
density. 



4.2 Log Priors 



In the case where we employed log priors, a far smaller number (and fraction) of models survive the 
confrontation with the experimental constraints. In this case, we randomly sampled 2 • 10^ models 
and found that only ~ 3000 models remained after imposing the theoretical and experimental 
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constraints. The main reason for this is rather straightforward to understand. In the log prior 
scan, although the range of the SUSY masses has been expanded, most of the generated parameter 
space points correspond to models with very light SUSY particles which are more highly constrained 
by existing experimental data and, consequently, far fewer of the models survive. Recall that we 
employed a smaller model sample in this case as our primary goal in undertaking the log prior 
study is to examine any differences with the flat prior analysis and consequently our presentation 
of the corresponding results will be somewhat brief. 

We begin our comparison with the flat prior results by examining the various Higgs and 
sparticle spectra. Figures [13] - [18] show the set of spectrum results for the set of log prior models 
which survive all of the constraints. These should be compared to Figs. [1]- [6] above in the case 
of flat priors. As can be seen these spectra are qualitatively very similar to what we have seen 
in the case of flat priors above, except for the obvious two effects already mentioned: the ranges 
of the allowed mass values have been extended so that larger masses are reached and the sample 
size is significantly smaller. However, the broad features discussed above for the flat prior sample 
also hold true for the log prior sample. The main feature that we observed in both log and flat 
prior samples is the apparent preference for relatively light spartners with this preference being 
somewhat stronger in the log prior case as would be expected. While the low end of the log prior 
distributions are similar to those in the flat prior analysis, naturally, the log prior spectra have 
longish tails which stretch out to significantly larger mass values. Larger A-terms, for example, can 
lead to greater stop mixing effects which then produce larger values for the light Higgs mass; this 
also leads to a slight increase in the maximum value of the Higgs mass in comparison to the flat 
prior case, as is observed in Fig. [THl 

Figures. [19] and [20] describe the nature of the nLSP and the corresponding mass splitting 
with the LSP for the case of log priors, as well as the associated electroweak origin of the LSP. While 
different in some details from the flat prior results, here we see that the properties of the LSP and 
the possible identities of the nLSP are essentially the same in the two cases. The slight difference 
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Figure 13: Distribution of slepton masses for the set of log prior models satisfying all of our 
constraints: selectrons/smuons(top left), staus(top right) and sneutrinos(bottom). 
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Figure 14: Distribution of squark masses for the set of log prior models satisfying all of our con- 
straints: 'light' (first /second generation) squarks are in the upper left panel, stops in the upper 
right panel and sbottoms are in the lower panel. 
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Figure 15: Distribution of neutralino and chargino masses for the set of log prior models satisfying 
all of our constraints. 
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Figure 16: Distribution of gluino masses and a comparison of the gluino and LSP masses for the 
set of log prior models satisfying all of our constraints. 
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Figure 17: Distribution of heavy and charged Higgs masses for the set of log prior models satisfying 
all of our constraints. 
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Figure 18: Distribution of the light Higgs mass on a log scale for the set of log prior models satisfying 
all of our constraints. 

in the identity of the nLSP between the log and flat prior scenarios is due to the difference in the 
statistics of the two sample sizes. Overall, the main distinguishing feature that we observe arises 
again due to the preference for lighter sparticle masses in the case of log priors and the reduced 
statistical power of the surviving sample. Complementary information on the Wino/Higgsino/Bino 
content of the LSP in the case of log priors is provided in Figure [2T1 Here we see results which are 
qualitatively similar to those found in the case of flat priors. 

The log prior predictions for a number of experimental observables are displayed in Figs. [22] 
and [23l These are seen to be quite similar qualitatively to those obtained in the case of flat priors. 
Some minor differences are: (i) the slight preference for a larger branching fraction for the 6 — > 57 
transition in the case of log priors; (ii) the lack of a bimodal peaking structure for A(g — 2)^ in the 
log prior case, which is directly linked to the lighter spartner masses and (iii) though peaked near 
the same value as in the flat prior sample, the relative number of models in the log prior sample 
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Figure 19: (Top) Identity of the nLSP. (Bottom) nLSP-LSP mass splitting as a function of the LSP 
mass. Both apply in the case of log priors. 
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;ure 20: Wino/Higgsino/Bino content of the LSP in the case of log priors. 
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Figure 21: Fractional Wino/Higgsino/bino content of the LSP in the case of log priors as functions 
of the LSP mass(top) and the nLSP-LSP mass splitting(bottom). Note that each of the log prior 
models has three entries on these figures. 
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Figure 22: Distributions of predictions for several observables as well as tan (3 for our model sample 
subject to the constraints discussed in the text in the case of log priors. The blue and green 
dashed lines show the SM predictions as well as the current central values obtained by experiment, 
respectively. 
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that predict larger values of Jl/i^ is found to be somewhat larger. All of these small effects can be 
traced to the log tilted mass spectra and the corresponding extended mass range. In Fig. [23] we 
see that the log prior models are perhaps even less fine tuned than those obtained in the flat prior 
case which is not a surprise since lighter spartner masses are now favored. 

From the analysis presented above, we conclude that on the whole it would seem safe to 
say that the flat and log priors lead to very similar results for almost all of the quantities examined 
with no significantly large differences. 

5 Discussion 

There are many ways to begin a detailed collider study of the varied signatures arising from the 
very large set of models that we have generated. Such studies are beyond the scope of the present 
paper but there are some observations that we can make in addition to those above. 

One of the best ways to probe the phenomenological implications of our set of models for 
the LHC is to examine the mass spectrum beyond the LSP and nLSP. Feldman, Liu and Nath 
(FLN)[65j have performed a previous study within mSUGRA where they categorize and examine 
the 22 possible orderings for the 4 lightest new particles in the spectrum. These are rank ordered by 
the number of times they occur while scanning the parameter space. This mass spectrum ordering 
is important as it influences the various cascade decay chains for heavy colored SUSY particles 
produced at the LHC. To get a feeling for our results, we have repeated this analysis using both 
our flat and log prior model sets. The first question we addressed was how frequently do we obtain 
a model which can be identified as one of the 22 found by FLN for either set of pMSSM priors. 
This is answered in Table [21 Interestingly, we see that one of the FLN orderings, mSPlO, does 
not occur in either of our flat or log prior model sample for the pMSSM at this level of statistics 
indicating that they occur very rarely. Many of the other FLN models are seen to occur rather 
rarely in our two model pMSSM sets, whereas the most common FLN models are also common in 
our two model sets. 
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Figure 23: Distributions of relevant dark matter predictions for our model sample in the case of 
log priors subject to the constraints discussed in the text. All cross sections are in pb. 
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Figure 24: The distribution of fine-tuning (= A) in our set of models with log priors and as a 
function of the LSP mass. Also shown is the correlation of the amount of fine-tuning with the 
predicted dark matter density. 
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% of Models 


mSP 


Mass Pattern 


Linear Priors 


Log Priors 


mSPl 


Xi<xt <X2< Xs 


9.82 


18.59 


mSP2 


X? < xf < < A/H 


2.08 


0.68 


mSP3 


Xi<Xi <X2< n 


5.31 


6.64 


mSP4 


Xi<xt <X2<9 


2.96 


3.73 


mSP5 


Xi<ri<iR< Dr 


0.02 


0.14 


mSP6 


Xi < n < xf < xl 


0.46 


1.22 


mSP7 


Xi<fi<iR< xf 


0.02 


0.03 


mSP8 


Xi <Tl < Ar^ H 


0.10 





mSP9 


X? < n < £r < A/H 


0.01 





mSPlO 


X? < fi < ii < Ir 








mSPll 


Xi < h <Xi <X2 


0.09 





mSP12 


n ~ -1- 

X\<ti<fi< xf 


0.01 





mSP13 


Xi < ii < fi < 4 


0.01 





mSP14 


X? < A ~ < 


0.35 


0.10 


mSP15 


xl<A^H<xf 


0.08 





mSP16 


~n ATT- 

Xi < Ar^ H <Tl 


0.01 


0.03 


mSP17 


Xi < n < X2 < xf 


0.18 


0.41 


mSPlS 


Xi<n<iR< h 


0.01 





mSP19 


Xi<Ti< h < xf 


0.01 





mSP20 


Xi < ii <X2<xf 


0.06 





mSP21 


X? < ii <fi<X2 


0.01 





mSP22 


Xi < X2 < xf < 5 


0.27 


0.51 



Table 2: Frequency of occurrence for the FLN|65j model sequences in our our flat and log prior 
pMSSM model samples. 
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The next question we addressed was how many mass patterns do we find for the 4 Hghtest 
new particles for both sets of priors and how common are they. This is a generaUzation of the FLN 
analysis to the pMSSM. In the case of flat (log) priors we find 1109(267) different mass patterns for 
the four lightest sparticles; the difference in these two numbers is largely due to the much greater 
statistics in the the flat prior case. However, in either scenario this number is many times larger 
than that found by FLN and it is likely that even more patterns may emerge if greater statistics 
were to become available. This result should not be surprising since we are imposing no constraints 
on the SUSY breaking mechanism and hence the sparticle mass ordering above the LSP is in general 
completely arbitrary. Table [3] shows the set of 25 particle mass orderings for the 4 lightest sparticles 
that occur most frequently amongst our surviving models for both flat and log priors. Note that 
the sets for the flat and log priors are generically similar but are different in detail. Note that the 
most frequently occurring FLN orderings are also among the 25 that we obtain for either set of 
priors. Conversely, however, most of the FLN orderings do not appear amongst our top 25. We 
note, however, that once the assumptions of minimal SUGRA are relaxed more of the new patterns 
that we find in the pMSSM begin to emerge. Four of these are given by the NUSPl-3 and DBSP6 
models as found in the last two papers of Ref. [65] . 

In conclusion, we see that the predictions of mSUGRA are greatly extended when examining 
the fuh pMSSM. 

6 Conclusions 

In this paper we begin a detailed study of the 19 parameter, CP-conserving pMSSM under the as- 
sumptions of (i) minimal flavor violation, (ii) the LSP being identified with the lightest neutralino, 
(Hi) only the Yukawa couplings of the third generation are important and {iv) the WMAP measure- 
ment of the dark matter relic density is taken only as an upper limit to that arising from the LSP. 
Carefully subjecting numerous combinations of parameter space points, with mass terms randomly 
generated assuming either flat or log priors, to a large number of theoretical and experimental 
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Linear Priors 


Log Priors 


Mass Pattern 


% of Models 


Mass Pattern 


% of Models 


y? < 


< y^ < y§ 


9.82 


X? 

A 1 


< 


xf 

A 1 


< y^] < y? 


18.59 


y? < 


< X9 < in 

^ A/ ^ ^it 


5.39 


X? 

A 1 
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xf 

A 1 


< Ya < Z>r 

A Z • 


7.72 


Y? < xt" 

A, 1 A, 1 


< v9 < fi 

^ A^ ^ ' 1 


5.31 


X? 
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xf 

A 1 


<X9 <iR 

• A Z ^ 
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X? 
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xf 

Al 


< Yo < Tl 

A 2 J- 


6.64 


X?<X^ 


< X9 < ^1 

Az J- 


4.89 


X? 

Al 


< 


xf 

Al 


< Xo <dR 

A Z -ft 


5.18 


X?<X^ 


< Xq < (ii? 

^ A.^ 


4.49 


X? 

Al 


< 


xf 

Al 


< Xo < i^f 


4.50 


x?<xf 

Al A-l 


< Xo <uri 

^ A.^ ^ "'JX 


3.82 


X? 

Al 


< 


xf 


< Yo < fei 

A z 1 


3.76 


Y? < 

/vl A. 1 


< Yo < 

A_^ C/ 


2.96 


X? 

Al 


< 


xf 


<X2<5 


3.73 


y? < 

A 1 A 1 


< v9 < 

^ Az ^ t 


2.67 


X? 

A 1 


< 


xf 


< X2 < ^i? 


2.74 


y? < y^' 

A 1 A 1 


< Ya < -Ur 

^ A^ 1j 


2.35 


X? 

Al 


< 


xf 


< i>T- < fi 


2.27 


X? < xf 


< Ur < Ti 


2.19 


X? 


< 


xl 


< xf < Xs 


2.24 


X?<X^ 


< xf < X3 


2.15 


X? 

Al 


< 


xf 


< 4 < X^ 


1.42 


Xi <xf 


<X2<^ 


2.00 


X? 


< 


xf 




1.32 


Xi <xf 
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1.40 


X? 


< 


Tl < xf < X2 


1.22 


Xi <xf 


<i>e<iL 


1.37 


Xi 

Al 


< 


xf 


< n < X2 


1.19 


Xi < xt 


<h<x2 


1.35 


X? 


< 


X2 < xf < i>r 


1.15 


Xi < xf 


< in < xl 


1.32 


X? 

A 1 


< 


iR 


< xf < X2 


1.05 


^ < < i7± < X? 


1.24 


X? 

A 1 


< 




< n < xf 


1.02 


X? < xf < < X2 


1.03 


X? 


< 


xf 


<i>e<h 


0.95 


Xi <xf 


<UL < di 


0.95 


X? 


< 


xf < 4 < X2 


0.71 


Xi < xf < ^1 < X2 


0.89 


X? 


< 




< xf < X2 


0.68 


x?<xf 


<UR<X2 


0.84 


X? 
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xf 


<xl<A 


0.64 


Xi <xf 


<A<H 


0.74 


X? 


< 


xf 


<^r<X2 


0.61 


Xi < xf 


<9<X2 


0.65 


X? 


< 


x^ 


<xf < dR 


0.54 


Xi < xf 


< fi < l>r 


0.51 


xS 


< 


xf 


<fi < Ut 


0.54 



Table 3: The 25 most common mass orderings for the four lightest sparticles amongst our pMSSM 
model sets assuming either flat (left) or log(right) priors. 
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constraints we arrive at two sets of models that survive all of the restrictions. 

With these two sets of models, we first examined the resulting distributions of expected 
Higgs and sparticle masses; qualitatively and even semi-quantitatively either choice of prior was 
found to return quite similar results except for trivially obvious differences due to how the param- 
eters were generated. In both cases, we found some instances where surprisingly light sparticles 
are still allowed by the data. Next, we studied the identity of the nLSP and its mass splitting with 
the LSP. We found that for either choice of priors, over a dozen sparticles can play the role of the 
nLSP with more or less the same probabilities in the two model samples. The most common nLSP 
identities were found to be xf and X2 with other possibilities being Cij, fi and ii. Less familiar 
nLSP identities were also found to occur such as the gluino or the first/second generation squarks. 
We also showed that the mass splitting between these two states, Am = rUnLSP — "mLSP, which 
is a critical parameter for collider signals and studies, was found to vary over a very large range 
of six orders of magnitude with values in excess of ~ 100 GeV and as small as ~ 100 keV or less. 
The greatest range was found in the case where the lightest chargino was the nLSP. Small values of 
Am were found to be quite common particularly in the case of the lightest chargino or the second 
lightest neutralino being the nLSP. This would imply that searches for long-lived and (essentially) 
stable particles at the LHC will be of particular importance. In some cases, squarks and gluinos 
were found to be relatively light and not widely separated in mass from the corresponding LSP. The 
decays of such states would likely lead to soft jets in the final state which may be difficult to observe 
above SM backgrounds at the LHC. The mass spectra of the models we obtained were such as to 
substantially diminish the needed amount of fine tuning in comparison to, e.g., mSUGRA-based 
scenarios. 

Subsequently, we examined the postdictions arising from both model sets for a number of 
observables such as the branching fraction for 6 ^ 57 and the shift in the value of the g — 2 of the 
muon; these were found to be reasonably similar for both cases of priors. Since the predicted dark 
matter density was allowed to float to any value at or below that obtained by WMAP, we found 
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that the average contribution of the LSP to dark matter peaked at ~ 4%, with long tails, implying 
that most of the time some other source of dark matter must typically dominate. However, we do 
see that there is a reasonable fraction of models where the LSP contribution to the dark matter 
density is, in fact, quite appreciable in comparison to the WMAP measurement. These results 
were found to hold for either choice of priors. We then went further and examined the implications 
of these models for the direct detection of dark matter; we found that the range of possible cross 
sections is far larger than that obtained in any of the well-known SUSY breaking frameworks, such 
as the cMSSM or mSUGRA. Furthermore these ranges extended into the regions one would have 
expected to only arise from either the Universal Extra Dimensions or Little Higgs models. This 
implies that these observables alone cannot be used to discriminate between different Terascale 
models. 

It is important to note that while the two sets of priors have led to qualitatively very 
similar pictures, some quantitative differences do appear as can be seen from the figures above. 
For example, the mass spectra of the SUSY particles in both cases, while qualitatively similar, are 
different in detail since the log prior case generally leads to smaller masses. 

Clearly the detailed study of the pMSSM models is quite exciting as well as complex; we 
are just at the initial stages of this analysis and there are many things remaining to be done which 
we will address in subsequent publications. In addition, we plan to make our model sets public in 
the near future, in order to encourage a wider set of possible analyses. Hopefully we will not have 
to wait too long until we actually have real data from the LHC to examine. 
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